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The azimuthal angular dependence of the nuclear modification factor Raa{pt , (p, Npa,rt) recently 
obtained by PHENIX is related at low pr to the trigger <j> dependence of the ridge yield as measured 
by STAR in a framework in which the azimuthal anisotropy is driven by semihard scattering near 
the surface. Careful consideration of the initial geometry leads to the determination of a surface 
segment in which the production of semihard partons are responsible for the <f) dependence of the 
inclusive distribution on the one hand, and for the angular correlation in ridge phenomenology 
on the other. With V2 also being well reproduced along with Raa and ridge yield, all relevant (f} 
dependencies in heavy-ion collisions can now be understood in a unified description that emphasizes 
the ridge production whether or not a trigger is used. 
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I. INTRODUCTION 

Recent measurement of tt" production in heavy-ion 
collisions, expressed in terms of the nuclear modifica- 
tion factor Raa, shows interesting dependence on the 
azimuthal angle 4> at various centralities and transverse 
momenta px [U]. Those properties in are lost when 
the azimuthal anisotropy is summarized in terms of the 
elliptic flow coefficient V2. Similarly, the production of 
ridge in association with triggers has revealing behavior 
in , the azimuthal angle of the trigger relative to reac- 
tion plane, at various centralities [2!|, but those properties 
are also lost upon integration over all (j)s in the determi- 
nation of the total ridge yield 0, |3]- In this paper we 
show that these two phenomena are related, even though 
RaaIpt, 4>, ^part) is a measure of the single-particle dis- 
tribution, while the ridge yield l^^(0s, ^part) is a mea- 
sure of the correlation between trigger and associated 
particles in specific pr ranges. 

In [l| the attention is given mainly to hard processes 
at high pt, but it is noted that at low px and large A0 
around 7r/2 (the (j) angle relative to the reaction plane) 
Raa is nearly constant in A^part ■ Such a curious property 
calls for an explanation. A difference in the tempera- 
tures along the x- and y-directions in the transverse plane 
has been discussed in the Buda-Lund model [1, Q where 
V2 can be expressed in terms of a scaling variable that 
depends on the difference between those temperatures. 
However, being a hydrodynamical study it does not con- 
sider explicitly the dependence. Azimuthal anisotropy 
at low Pt that does not depend on the assumption of 
rapid thermalization has been investigated in the con- 
text of semihard scattering and the ridges that are pro- 
duced 0, [1] ■ There also V2 is derived without detailed 
consideration of the (p variable. The azimuthal correla- 
tion between the ridge and trigger particles has, however, 
been studied in the correlated emission model (CEM) Q . 
It is from that study that we shall find in this paper an 



explanation of the cf) dependence on Raa- 

Since Raa was introduced originally for the purpose of 
exhibiting the effect of jet quenching at high pT, it is ac- 
tually inappropriate to describe the behavior at low px- 
In the next section we shall introduce a modified Raa 
that is more suitable. The present paper is focused ex- 
clusively on Pt < 2 GeV/c. It is conventionally thought 
that at such low pT only soft physics is involved, a view 
that has been regarded as being too simplistic on account 
of the semihard scattering that can enhance the soft ther- 
mal partons [l^, [HI . Whether or not trigger particles 
are used to select a subset of events, semihard partons 
are pervasive and generate ridges that affect both single- 
particle and dihadron distributions. That is why we are 
able to relate the two azimuthal features of the hadronic 
observables that have been determined by the two RHIC 
experiments [U, 

The focus of this paper is on the azimuthal de- 
pendencies of RAA{PT,4>,N.p.^^t) and the ridge yield 
Y^{(j)s,Npa.Tt) at midrapidity, where data exist [U, 0]. 
We leave out completely any consideration about lon- 
gitudinal correlation, on which there is experimental ev- 
idence that the ridge particles assoicated with a trigger 
can be widely separated in rapidity with A77 ^ 4 [12| . 
although autocorrelation without trigger shows a ridge 
with A77 PS 2 [1^. The physics for longitudinal correla- 
tion is very different from that of azimuthal correlation 
at midrapidity. We regard them as separate issues to be 
treated separately, though all basic problems are related 
at some deep level that we do not address here. 



II. CONNECTION BETWEEN Raa (prA^Np) 
AND THE RIDGE 

The main content of Raa{pt:<I>tNp) is the single- 
particle distribution of AA collisions, which we abbre- 
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viate as pi: 



Pt dpT d(p 



(1) 



where Np is a shortened notation for the number of par- 
ticipants, A'^part- In the foHowing we shall consider only 
the production of tt", since that is what is measured in [1] 
with high accuracy. The properties of pi{pT, 4>, Np) are 
usually presented, especially in hydrodynamical studies, 
in terms of the px spectra (average of pi over 0) and 
the elliptical flow coefficient V2{pt) (second harmonic in 
(f>) . The conceptual basis of flow deemphasizes the role of 
hard scattering among partons and offers a fluid descrip- 
tion of the bulk medium at low px- In that description 
azimuthal anisotropy in non-central collisions is regarded 
as the consequence of asymmetric pressure gradient in the 
transverse plane at early time, assuming fast thermaliza- 
tion. However, that is not the only approach that can 
claim physical relevance, although it is the conventional 
view backed by a large number of investigations based on 
hydrodynamics [T^ . 

An alternative view is to regard semihard scattering 
near the surface of the nuclear overlap as the driving force 
of the azimuthal anisotropy, since at low enough virtu- 
ality such scattering processes are pervasive throughout 
the medium, and when they occur near the surface of 
the almond-shaped initial configuration of the dense sys- 
tem, the semihard partons can emerge from the medium 
and not only hadronize as intermediate-pr jets, but also 
generate ridge particles at lower pT [ill- Many is- 
sues are involved in the above statement. First, the 
fact that semihard scattering cannot be calculated re- 
liably in pQCD does not mean that its effects are not 
important. Second, those effects are sensitive to the 
geometrical shape of the initial system at early time, 
since if the parton's transverse momentum is kx ~ 2- 
3 GeV/c, the time scale involved is ^ O.I fm/c. The 
semihard jets created near the surface give rise to non- 
trivial dependence of low-p-r partonic distribution in- 
dependent of the validity of the notion of pressure gra- 
dient at T < 1 fm/c 0. Third, those ^-dependent soft 
partons hadronize and form the ridge structure observed 
by STAR whether with trigger llBl or without trig- 
ger |l3|. A theoretical description of the connection be- 
tween ridges and V2 in the approximation of using a very 
simple geometrical picture is given in 0, Q . It is along 
the same line of reasoning, but using a more realistic 
treatment of the initial geometry, that we relate in this 
paper all features of 0-dependent observables, namely: 
Raa{pt,(I>,Np), V2{pT,Np) and ridge yield Y'^{4>s,Np) 
for PT < 2 GeV/c. 

For hard scattering calculable in pQCD it has been 
shown that the Landau-Pomerachuk-Migdal interference 
effect suppresses the radiative energy loss of a hard par- 
ton in the initial phase of its trajectory [H, [l3|- Apply- 
ing the inference to semihard partons, despite uncertain 
validity, there can be an initial time interval when no en- 
ergy loss may occur, the consequence of which is that the 



creation points of the relevant semihard partons can be 
farther away from the surface. What count for the ridge 
formation are those semihard partons that do lose en- 
ergy on the way out and drive the azimuthal anisotropy. 
A shift in the time period when that occurs in the tra- 
jectory of a parton in the medium does not alter the re- 
lationship we propose between semihard scattering and 
ridge formation. 

An important point to stress in distinguishing our ap- 
proach from the conventional hydro approach is that the 
bulk of what is calculated in the latter contains <f> depen- 
dence at Pt < 2 GeV/c and serves as the background of 
higher pT hard processes, while the bulk in our approach 
has no cf) dependence and it is the ridge component at 
Pt < "2 GeV/c that gives rise to the 4> dependence of the 
soft component. Thus we write 



pi(pT, (t>, Np) = B{pT,Np) + R{pT, (f>, Np) 



(2) 



with B and R denoting bulk and ridge, respectively, for 
PT < 2 GeV/c. 

At a fixed pt, Pi(pt, 4>7 Np) increases rapidly with Np, 
but Raa{pt, 4>, Np) decreases with Np except when (j) 
7r/2. The relationship between them is 

RAA{PT,(t),Np) = Z{pT,<l>,Np)pi{pT,(l),Np), (3) 

where the rescaling factor is 



Z(pT,0,iVp) = [N^{Np)dNpp/pTdpj 



(4) 



Nc being the number of binary collisions, short for A'coii, 
and dNpp/pTdpTd(j) the single tt'^ inclusive distribution 
in pp collisions. This rescaling factor is designed to give 
Raa a quantitative description of the effect of jet quench- 
ing at high Pt. But at low pT, which is the region of 
interest to us in this paper, the more relevant rescaling 
factor is 



Z'{pT,(j),Np) = [NpdNpp/prdpq 



(5) 



where Np is used instead of Nc because of the dominance 
of soft processes. Thus we define a different nuclear mod- 
ification factor 

R'aa{pt, 0, Np) = Z'ipT, Np)pi{pT, 0, Np). (6) 

From the data on Raa given in P], we can determine 
^AA- We illustrate their differences by showing in Fig. 
1 (a) the original Raa(pt,4'i Np) and (b) the modified 
R'^j^{pt,4': Np), for (j) in the ranges < Acj) < 15° and 
75° < A(j)< 90° for I.O < pt < 1.5 GeV/c. In the follow- 
ing we shall for brevity use cj) = and 7r/2 to denote the 
two angular ranges, although in calculations that involve 
comparison with data we shall use (p — 7.5° and 82.5°, re- 
spectively. The error bars in Fig. 1(b) are calculated from 
the errors on Np and N^ neglecting the errors on Raa 
which are much smaller. Comparison can also be made 
for 1.5 < Pt < 2.0 GeV/c, which we omit here since they 
are similar to those in Fig. 1. The effective values of the 
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two pt ranges will be 1.21 and 1.71 GeV/c given numer- 
ically in tables in We see in Fig. 1(a) that Raa for 
(j) = decreases rapidly with Np, but R'^a (b) is nearly 
independent of Np for (p — 0. Whereas Raa at (j) — Tr/2 
is nearly constant, — Tr/2) increases linearly with 

Np. Intuitively, these properties of the Np dependence for 
R'aa ^^'^ easier to understand, since the elliptical shape 
of the initial nuclear overlap is relatively independent of 
the impact parameter b when viewed from the broad side 
at (j) = 0, but becomes narrower as b increases, when 
viewed from = 7r/2 where one sees mainly the narrow 
side of the ellipse. Our problem is to reproduce these 
properties of R'aa quantitatively in a framework where 
the decomposition of pi into the two components in Eq. 
^ plays an important role. 

As described in [7, 8] and summarized in the beginning 
of this section, the azimuthal anisotropy arises from co- 
pious semihard scattering near the overlap surface; the 
produced semihard partons lose energy to the medium, 
whose enhanced thermal partons hadronize into the ridge 
particles that carry the footprint of the initial geometry. 
The (j) dependence of the ridge distribution R{pt, 0, Np) 
is therefore expected to be minimal at large Np (small 
&), but increases with decreasing Np, thus opening the 
gap between cj) — and 7r/2 as shown in Fig. 1. The re- 
coil semihard partons directed toward the interior of the 
dense medium are thermalized and retain no memory of 
the initial geometry by the time they are hadronized; 
together with all the other soft processes that occur in 
the interior they form the bulk that is independent. 
With B{pt, Np) in Eq. ^ being independent of (f>, it is 
possible to calculate V2{pT,Np), whose value depends on 
the relative magnitudes of B{pT,Np) and R{pT,4>,Np). 
The absolute magnitude of R{pt, 0, Np) can be related 
to the ridge yield that depends on the trigger angle (ftg in 
dihadron correlation 0]. Thus we have a tightly con- 
strained system that restricts the options available to 
finding a satisfactory physical basis to explain all the 
relevant experimental features. 

To summarize this section, we have described an over- 
all picture in which the (f> and Np dependencies of 
Raa{pt,4':Np) are to be related to those of the ridge 
structure at low pt- In the following sections we imple- 
ment this picture with concrete calculations. 



III. RIDGE AND DIHADRON CORRELATION 

The discussion in the preceding section is concerned 
mainly with the single- particle distribution pi{pT, (f>, Np). 
However, ridgeology (short for ridge phenomenology) in- 
volves two-particle correlation. Let us denote two-pion 
distribution after background subtraction by 



P2(Pl,01,P2,02) 



dNll 



Pi^dpi^d(f>ip2j,dp2Td(l}2 



(N,), (7) 



where the subscript T will for brevity be omitted, since 
all momenta will hereafter be restricted to the transverse 
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FIG. 1: (a) Original Raa and (b) modified R'aa ^t fixed pr in 
the range 1.0 < pr < 1-5 GeV/c for (j> in the ranges < A<f) < 
15°, denoted hy (/> ^ (in solid circle), and 75° < A4> < 90°, 
denoted by = 7r/2 (in open squares). The data are from [J. 



plane with \ri\ < 0.35. We shall use label 1 to refer to 
the trigger and label 2 to a particle in the ridge. Thus, 
we have A(j> = (t>2 — <t)i for the correlation angle, not to 
be confused with the notation in [l| where Acj) denotes 
the angle of a pion relative to the reaction plane. When 
the reaction plane is used as the reference for triggers 
in correlations (as will be done in the following), then 
4>i is identified with (ps in 0]. The dihadron correlation 
distribution has two components, which are referred to 
as Jet and Ridge 0, 0] , the former being the peak in A77 
that sits above the pedestal that is the latter with a wide 
spread in A77. We write them without the arguments as 



P2 = pi + P2- 



(8) 



We shall consider only the latter, assuming that in rid- 
geology the former has been subtracted out from the di- 
hadron correlation function. 

The stage is now set for us to discuss the crucial point 
in our physical basis that relates p^{pi., 0i,P2, ^2) to the 
ridge component R{pt, 0, Np) in Eq. ^ for the single- 
particle distribution. Usual ridgeology involves the use 
of a trigger and measures the momentum of another par- 
ticle associated with it. For pi in the range 3-4 GeV/c 
as studied by STAR in Q , the parton that produces the 
trigger is semihard, i.e., its transverse momentum is in 
the intermediate kx region, < 5 GeV/c. PHENIX uses 
even lower trigger momentum (2 < pi < 3 GeV/c) to 
study the away-side structure [l8|. The azimuthal prop- 
erties of the ridge structure as functions of 0i and (f>2 
separately provide significant insight into the correlation 
between the semihard and the soft thermal partons that 
does not depend sensitively on the trigger momentum Q . 
A model on that azimuthal correlation has successfully 
reproduced the data on the (j)s dependence of the ridge 
yield Q; moreover, a prediction on the asymmetry of the 
ridge structure in that model (GEM) has recently been 
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verified by an analysis of tlie STAR data . Tlie signif- 
icance of tliat finding as related to our present problem is 
that the ridge yield studied in 0, Q involves an explicit 
integration of the Ai/) distribution of over the ridge 
angle 4>2, while the ridge component R{pt, 4>, Np) in pi 
involves an implicit integration of the same over the trig- 
ger angle 4>i. A successful description of the correlation 
between (pi and 02 is therefore a crucial input for the 
determination of R{pt, Np). 

An important part of the physics implied in the above 
is that the direction of the semihard parton that gener- 
ates the ridge need not be identified by a trigger particle, 
i.e., a ridge can exist whether or not there is a trigger. 
A ridge particle at (j)2 can be due to semihard partons at 
a range of 0i around (f)2, so for an inclusive ridge distri- 
bution in (f)2, the value of 4>i should be integrated over, 
relative to every fixed 4>2- If we omit the arguments px 
and Np in R{pt, 4>, Np) and identify 4> = 4>2, then we have 

i?(02)cx J d0i pf(01,02), (9) 

where the proportionality factor involves the integration 
over pi, which will not be considered explicitly. On the 
other hand, the ridge yield determined in Q integrates 
over all 02 in the range |A0| < 1, so we have 

r«(0l)cx J d02 pf(01,02), (10) 

where the proportionality depends on the number of trig- 
gered events in prescribed ranges of pi and p2, since 
is the per-trigger yield. The effective range of integra- 
tion in both Eqs. ^ and (fTO)) depends on the correlation 
width in p|^(0i, 02) and is much narrower than 1 rad. 

It should now be clear why the dependence of Raa 
in the PHENIX data [l| should in our view be related to 
the ridgeological study by STAR 0. There is, however, 
more complication than just the correlation in A0. The 
dependence on ccntrality that has not been made explicit 
in Eqs. (0) and ^ is evident in both Q and [2]. To 
achieve the correct dependence on Np requires a careful 
consideration of the geometry of the problem, which is 
the subject of the next section. 

A summary of this section is that the distribution of 
ridge particles and the ridge yield are two sides of the 
same coin, the former being the projection onto 02 of the 
two-particle distribution /9|*'(0i, 02), the latter being the 
projection of the same onto 0i. 

IV. GEOMETRICAL CONSIDERATIONS 

The correlation between a semihard parton at 0i and 
a ridge particle at 02 has been studied in the GEM that 
involves several subprocesses (a) The semihard par- 
ton loses energy to the medium, (b) Successive radiation 
by the parton enhances the thermal partons in the vicin- 
ity of the trajectory, (c) The enhanced thermal partons 



are carried by local flow to various points on the surface 
in directions normal to the surface, (d) Hadronization of 
those enhanced thermal partons by recombination gen- 
erates the ridge particles that are centered around the 
average flow direction, (e) Density of the ridge is highest 
when the fiow directions are in alignment with the direc- 
tion of the semihard parton. (f ) The resultant correlation 
between 0i and 02 can be described by the convolution 
of a Gaussian distribution in 0i — ?/> (where "0 is the lo- 
cal flow direction) with another distribution describing 
the fluctuation of 02 from ip. In the midst of all those 
subprocesses that cannot be calculated from first princi- 
ples, there are considerations of nuclear density, points of 
semihard scattering, emission points along the trajecto- 
ries, and the exit points at the surface, all of which have 
been taken into account to the extent possible. What we 
can adapt from that study in ridgeology is the central 
result that there is a correlation 

r<t A \ J -'^ix,y)f\ 
C(a;,y,0i) = exp ^ > (11) 

between the semihard parton direction 0i and the flow 
direction ip{x, y) at the exit point (x, y) on the surface, 
with (7 = 0.33. The angle between the ridge particle 
02 and ip due to fluctuations has negligible effect on the 
results. 

For every semihard parton direction 0i a variation in 
the starting point of the trajectory in the nuclear overlap 
leads to a variation in ip{x,y). This important property 
can be turned around and applied to the ridge component 
of the single-particle distribution. With 02 identified as 
■0, for every ridge particle detected at 02 there is a variety 
of semihard partons in direction 0i that can contribute 
to that ridge depending on where the exit point (x, y) is, 
so long as 101 — 02 1 is within the Gaussian width of Eq. 
PT|) . The point to stress here is that for R{px,4>,Np) 
in the single-particle distribution in Eq. is 02 and 
the semihard parton at 0i is undetected and should be 
integrated over. 

The scheme described above can be implemented by 
concrete calculation because the main input is angular 
correlation using early time geometry. Hadronization, 
of course, occurs at late time, and we do not rely on 
hydrodynamics to describe the evolution of the system 
in space-time. Our assumption is that the ridge particles 
in the final state have the same distribution as that 
of the enhanced thermal partons at early time, which in 
turn is prescribed by Eq. (fTTI) . The bulk component that 
has no preferred direction to expand (not being the same 
as the hydrodynamical bulk) is isotopic, and becomes the 
background of the ridge that carries all the information 
about the anisotropy of the initial geometry. 

When two nuclei of radius Ra collide at impact param- 
eter b, the almond-shaped overlap has width and height 
given, respectively, by 

w = l^b/2, /i = (1 - 62/4)1/2, (12) 
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where all lengths are normalized by Ra ■ The appropriate 
geometry that can best describe the intial configuration 
almost immediately after impact is the ellipse 



w h 



(13) 



The flow direction at the surface is normal to that sur- 
face, so its azimuthal angle is 



V'(a;, y) — tan 



uP"y 



(14) 



It should be noted that although we refer to ip as the flow 
direction, the use of those words is mainly a compact way 
to name the direction normal to the surface, as defined in 
Eq. We place no emphasis on hydrodynamics as the 
basis for that description, especially in the initial phase 
of the evolution of the medium. The implication of Eq. 
PT|) is that for a fixed semihard parton at the ridge 
particle at (/)2 identified with ipi^, y) can deviate from (pi 
at most by cr, as illustrated in Fig. 2(a). By the same 
reasoning, for a fixed 4>2 the possible semihard partons 
that can contribute to the corresponding ridge must have 
their (pi that can deviate from (p2 by no more than a, as 
illustrated in Fig. 2(b). That means the exit point {x,y) 
of the semihard parton at the surface is restricted to a 
certain range, which we shall denote by S{(p2, b). 





(a) 



(b) 



FIG. 2; (Color online) Sketches to illustrate the relationship 
between the semihard parton's direction (j>i (in thin red ar- 
rows) and the ridge particle's direction p2 (in thick green 
arrows), (a) Three possible semihard partons all having the 
same pi, but generating ridge particles in different directions 
p2 not deviating by more than a from pi. (b) For a fixed 
p2 (two being illustrated), pi of contributing partons are re- 
stricted to within a cone of width a around p2 ■ 



From the geometry of the ellipse we can use the con- 
straint of Eq. (IT51) to write x = w cos and y = h sin 6 so 



that S{p2, b) can be determined by 

S{P,b) ^ f di^J [{dxf + {dy)^f' 
[w^ sin^ 9 + cos^ 9] ^''^ d0 



= h[E{02,a) - E{9i,a)] , (15) 
where E{'0, a) is the elliptic integral of the second kind 



Ei^,a)^ / (1 
Jo 



a sm 



)'/^d0. 



(16) 



In Eq. p5)) we have 
a{b) = 1 -2/^2 
6i = tan 



-1 



h 

— tan ( 

w 



i = 1,2, if (P, < tt/2 , 



If (p2 > tt/2, then 



cot 



— tan (p2 

w 



(17) 



(18) 



Thus the physical meaning of S{(p, b) is that it is the 
segment of the surface through which the semihard par- 
ton can be emitted to contribute to a ridge particle at 
p. In Fig. 3 we show S{p,b) for = and tt/2, and 
a — 0.33. The two curves start out with the same value 
at 6 = 0, but develop a wide gap between them, as b 
increases. If they are plotted against Np, then the be- 
havior of S{p, Np) appear similar to the Np dependence 
of R'jiA shown in Fig. 1(b) for b < 1.5, corresponding to 
Np > 50, apart from a common background for the two 
sets of data points. That is a strong hint that we are on 
the right track. 




FIG. 3: Surface segment S{p, b) vs normalized impact param- 
eter 6 for p = (solid line) and p = n/2 (dashed line) and 
a = 0.33. 

Note that the initial slope of S{0, &) at 6 = is pos- 
itive, but that of S{'K/2,b) is negative. The latter can 
intuitively be understood because the ellipse is narrower 
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than the circle when viewed from the top, so the segment 
of the surface within the width a decreases with increas- 
ing b. In the former case the elhpse is flatter than the 
circle when viewed from = 0, so at fixed a the seg- 
ment covers more of the surface as b is increased near 
b — 0. What is remarkable is that such detailed geomet- 
rical properties seem to be exhibited by the behavior of 
R'aa high Np. 

It should be mentioned that our consideration so far 
has been exclusively in the transverse plane at midrapid- 
ity so that the use of elliptic geometry has been simple 
and adequate for the calculation of the surface segment 
S{(f>, b) through which semihard partons are emitted at 
small \y\. Since semihard scattering occurs at early time 
(r < 0.2 fm/c), the appropriate geometry is the overlap 
in the transverse plane of two Lorentz-contracted discs. 
Longitudinal geometry enters the problem in the consid- 
eration of local density later in this section. 

We now assert as the basic assumption about the 
azimuthal anisotropy of the single particle distribution 
pi{pT,(t>,b) that it arises from the (j) dependence of the 
ridge component that is proportional to S{(j),b). That 
component depends on all three variables px, 4> and b. 
Since the thermal partons are exponential in px , we write 
in the factorizable form 



(19) 



where the use of the transverse kinetic energy Et is dis- 
cussed in as a means to account for the hadronic mass 
mo with Et{pt) = (Pt + "^o)^^^ ~ ^o- The factor A{b) 



depends only on b, since the other two factors in Eq. p9)) 
exhaust the pt and </> dependencies. The bulk term in 
Eq. © has no (f> dependence, so we write as in Q 



B{pT,b) 



-Et{pt)/T 



(20) 



where the prefactor /6 follows from the thermal parton 
distribution q'^dN^ / dqTdcp — Cqt exp{—qT/T) through 
recombination. C{b) is known from a previous study and 
is described as a power law in Np [20j , on which we shall 
make a slight modification below to render a better de- 
scription of the data. Since the ridge particles are due 
to the hadronization of the enhanced thermal partons, its 
T' is higher than the bulk T by a small amount discussed 
previously 0, Q and will be detailed below. 

The factorized form for R{pT,(j),b) in Eq. (|T9l) is an 
approximation of a complicated expression that involves 
an integration over all points of semihard scattering in 
the nuclear overlap region and another integration along 
the trajectory of each semihard parton in the medium 
0, as well as over the angle (pi in Eq. ([9]). The prob- 
abilities of the creation of a semihard parton and of its 
energy loss depend on the local density of the medium. 
We simplify all that by writing A{b) in Eq. (|19l) as being 
proportional to the effective density D{x' ,y' ,b) at a rep- 
resentative point (x' , y') near the surface where semihard 
scattering takes place. Factoring out C^(6)/6 explicitly 



which carries the dimension (GeV) ^ for both B{pT,b) 
and R{pt, (pi b), we write 



A{b)^ 



C\b) 
6 



aD{x',y',b) 



(21) 



where a is a free parameter to be fixed by the magni- 
tude of V2{pt, b) that determines the degree of azimuthal 
anisotropy. As we have noted earlier, this is a highly 
constrained problem in which many pieces of data are 
brought to bear on the phenomenology. 

For the local density D{x' , y', b) in the transverse plane 
we use the Glauber model to compute it, which, for an ar- 
bitrary point {x,y)^ is, apart from an unspecified overall 
normalization. 



D{x,y,b) = LA{x,y) 



1 - e 



-uLB{x,y) 



+LBix,y) [i _ £-'^^■^(^•2') 



(22) 



where uj — 4.6 for Au+Au collision ^] and La,b{x, y) are 
the longitudinal lengths 

LA,B{x,y) = — / dzp{s,z) , (23) 

PO J~ZA.B 

and p{r) is the nuclear density in Woods-Saxon form 



Pir) = Pq 



1 + e('—o)/« 



(24) 



with ro = 0.92, £, = 0.08, and po = 0.285 where all 
lengths are in units of Ra = 7 fm and p is normalized by 



A. In Eq. 



„2 



1 



1 — Is — 6r, and 



{x + b/2y + y 



(25) 



The approximation we make to avoid averaging over all 
points {x,y) in the initial overlap is to choose {x',y') to 
be at a representative distance Sr from the surface on the 
short axis of the almond. Detailed study has shown that 
the layer in the medium just below the surface in which 
semihard partons are predominantly produced, leading 
to the ridge particles, is of thickness 0.2 0. Thus we 
adopt Sr = 0.17 for the average and choose 



c' = 1 - 6/2 - Sr 



y'^0. 



(26) 



In this way D{x\y',b) is a calculable quantity for ev- 
ery b. Its normalization is absorbed into the unknown 
parameter a in Eq. (I2ip . 

The treatment of the geometry of the collision prob- 
lem in this section is far more realistic than the 0(0) 
function used in [tI, |1] . The incorporation of the corre- 
lation properties found in GEM '9'| in the determination 
of S{(j),b) endows the ridge component R{pT,4>,b) with 
a reliably calculable (j) dependence that can be tested as 
the physical origin of the anisotropy of the single-particle 
distribution pi{pT, </>, b). 

A way to summarize what has been done in this sec- 
tion is to use Fig. 2(b) and state that, for every ridge 
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particle at (f>2 indicated by the thick green arrow, the 
possible angles (pi of the semihard partons that can con- 
tribute to it are within the cone shown by the thin red ar- 
rows, which subtends a segment on the elliptical bound- 
ary. That segment is quantified by S{(j), b) that we have 
calculated. It gives the (p dependence of the ridge distri- 
bution R{pt, 0, whose h dependence is affected by the 
medium density that is also calculable from geometry. 



V. ELLIPTIC FLOW COEFFICIENT v2{pT,b) 

It is now straightforward to calculate U2, defined by 
the second harmonic of pi 



V2{pT,b) = (cos 2.; 



J^'' dcf) cos 2(t)pi {pT, <j), b) 



lo'' d.(j)pi{pT,(l), b) 
Using Eqs. ^ and ([TO)) we obtain 

j^'"^ d<j> cos 2(PS {(f), b) 



(27) 



V2ipT,b) = 



K{pT,b) + j;'^dct>s{<p,b) 



(28) 



where, owing to Eqs. pOt and ([2T 



K(pT,b) 



TrB{pT,b) 



2e-ST(PT)/T'^(^) 



2aDix',y',b) 



^Et(pt)/T" 



with 



1 



1 1 



(29) 



(30) 



Note that C^(6) in A{b) and B{pt, b) are cancelled, leav- 
ing only the density D{x', y', b) to prescribe the b depen- 
dence in K{px, b). 

The integrals in Eq. (|28l) are easy to evaluate, but the 
result on V2 depends on the magnitude of K{pT, b), which 
has an undetermined parameter a. Moreover, the suit- 
able values for the inverse slopes T and T' require some 
discussion of the thermal distributions. In the recombi- 



nation model for heavy-ion collisions [21[ the shower par- 
ton distribution S can be calculated from semihard and 
hard scattering, but the thermal parton distribution T is 
determined by phenomenology. It is the hadron distribu- 
tion at low pt that is fitted by an exponential form, from 
which is inferred the T distribution through TT recom- 
bination in the case of meson. TS recombination is then 
calculated to obtain the effect of semihard scattering at 
intermediate pr ■ What is done in is to extract the in- 
verse slope from the inclusive spectrum averaged over all 
(/). Ridgeology in 0, Q separates T' from T, but is based 
on a very simple model on (j) dependence. In our more 
realistic study here the complication of R{pt, b) in Eq. 
(IT^ cannot be isolated from the background B{pT, b) in 
Eq. ^ without the knowledge of the parameter a in Eq. 



(1211) , which in turn cannot be determined separately from 
T and T' . An iterative process may be necessary, the first 
step of which is what we do now by taking T' to be as 
given in [2l|: T' = 0.317 GeV, and AT = T' - T = 0.045 
GeV as given in Q. From Eq. §0^ follows T" = 1.916 
GeV. With these values we can calculate K{pT, b) and 
then V2{pt, b) in Eq. (P5|) with a suitable choice of a. 



0.3 



^ 0.2 



0.1 



PHENIX 



p^=1 .21 GeV/c 
p^=1 .71 GeV/c 




0.5 



1 

b 



1.5 



FIG. 4: Elliptic flow coeflicient V2{pT,b) vs b for pr = 1.21 
GeV/c (solid) and 1.71 GeV/c (dashed). Data are from [l[. 

In Fig. 4 we show our calculated results for V2{pT,b) 
for PT = 1.21 (solid line) and 1.71 GeV/c (dashed line) 
exhibiting excellent agreement with the data that have 
the corresponding (px) [J- The one parameter we adjust 
to get the best fit is a = 0.47. The nature of the b depen- 
dence is critically dependent on the properties of S{4), b), 
while K{pT,b) plays a less sensitive role, being of order 
1 for most values of & at pt < 2 GeV/c. K{pT,b) is a 
measure of the relative strength of the bulk component 
B{pT,b) to the ridge component R{pT,(l),b) in Eq. ^ 
for single-particle distribution. The fact that we have a 
good reproduction of the data on V2{pt, b) validates our 
approach to assigning the (f> dependence of pi(pT,0, &) 
entirely to R{pT, 4>,b). That validation is strengthened 
when the same formalism can be shown to reproduce the 
data on the ridge yield in 0s-dependent dihardron corre- 
lation. 

To summarize, we have in this section made the first 
phenomenological demonstration that the (j) dependence 
calculated in Sec. IV leads to a good description of the 
data on V2{pt, b) using one free parameter a for the over- 
all magnitude of the azimuthal anisotropy. 



VI. DEPENDENCE OF RIDGE YIELD ON 
TRIGGER DIRECTION 

We now consider the correlation between a trigger and 
the ridge particle associated with it. It may be helpful 
to recapitulate what we have done. For single-particle 
distribution the ridge component R{pt, 0, b) is related to 
the dihadron correlation distribution by Eq. ([SJ, where 
P2 is the ridge component of p2 expressed in Eq. ([5]). 
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Dominance of by the correlation function C{(j)i,(j)2) 
shown in Eq. ([TT|) with (j)2 = tp{x,y) leads to the study 
of the surface factor S'(0, 6) given by Eq. (fTSl) . The in- 
tegral in Eq. ^ is carried out in Eq. P3|) that results 
in R{(j)2,b) cx S{(f>2,b), as shown explicitly in Eq. ([T9|. 
Now, with the same P2^(0i, 02) we can calculate the ridge 
yield as given in Eq. PU)) . where the integration 

of P2{4'ij 'P2), dominated by C{4>i, 4)2), over 02 results in 



i,b) cx 5(01,6) 



(31) 



Thus what we know already about 5(0, 6), shown in Fig. 
3, need only be replotted to be compared to data. 

In the per-trigger ridge yield is analyzed for p^'^ 
and p^^°'^ in the narrow ranges 3 < p^'^ < 4 GeV/c 
and 1.5 < < 2.0 GeV/c in 6 bins of 0^ (trigger 

direction relative to the reaction plane) between and 
tt/2 for two centralities 0-5 % and 20-60 %. For the 
corresponding b values we use b — 0.33 for the former and 
average over 1.06, 1.24, 1.41, and 1.57 for the latter. From 
5(0, &) in Eq. ([T5|) we obtain the results for Y^{<l),b), 
shown by the solid lines in Fig. 5(a) and (b), compared 
to the data from [2] . The normalization of the theoretical 
curves is adjusted to render a good overall fit in Fig. 5(a), 
since the data on yield depend on the experimental cuts 
so their magnitudes are not predictable. However, once 
the normalization factor for Fig. 5(a) is determined, that 
for Fig. 5(b) is not readjustable. Our results show good 
agreement with the data both in the two 0s dependencies 
and in the normalization for the non-central collisions. 
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FIG. 5: Ridge yield per trigger vs trigger angle 0s relative to 
reaction plane for (a) 0-5% and (b) 20-60% centralities. Solid 
lines are calculated results from Eq. pip , and the data are 
from 0]. 

The mild decrease of Y^{(j)s) with 0s in Fig. 5(a) can 
be seen in Fig. 3 when b is restricted to 6 = 0.33. The 
precipitous decrease in Fig. 5(b) corresponds to the large 
gap between the two curves of 5(0, b) for = and 7r/2 
when 6 is between 1.06 and 1.57. Of particular interest is 
the increase of 5(0 = 0, b) when b is increased from 0.33 
to about 1.5. Indeed, the data of Y^{(j)s, b) in Fig. 5 show 



that the yield at the lowest 0s is larger for 20-60% than 
for 0-5%, unlike the situation for all higher 0s bins. That 
is a remarkable confirmation of the validity of Eq. ([3T|) 
and of the properties of 5(0, b) that we have obtained. 

Thus in this section we have shown the other side of 
the coin in the metaphor used at the end of Sec. III. The 
ridge yield is also well described by 5(0, b). 



VII. NUCLEAR MODIFICATION FACTOR 

Having shown the existence of substantial support for 
the relevance of 5(0, &) to V2{pT,b) and Y^{(j)s,b), we 
now return to RAAiPT-, 0, b), the new data on which stim- 
ulated our present line of investigation in the first place. 
Recall that in Sec. 2 an alternate nuclear modification 
factor (NMF) i?^^(pT, 0, Np) is defined by rescaling with 
Np instead of N^, for pt < 2 GeV/c; its dependence on Np 
is shown in Fig. 1(b). The way in which the data points 
for = and 7r/2 converge as Np approaches 400 is very 
similar to the way that 5(0, b) and 5(7r/2, b) converge as 
6 ^ in Fig. 3. The 0-dependent part is, however, only 
the ridge component in pi(pT,4>i^p)] the bulk compo- 
nent that is roughly of the same magnitude is the other 
piece yet to be added, rendering i?^^(pT, 0, ^p) to be- 
have as in Fig. 1. 

The normalizations of both B(pT,Np) and 
R{pT,<p,Np) in Eq. ([2]) should now be determined. 
Their relative magnitude has already been examined in 
connection with K{pT,b) in Eq. ([^ . Using Eqs. pT 
(PU)) and (HH), let us re-express their sum as follows 

piiPTA,Np) = C^iNp) fipT,^,Np) , 

where 

f{pT, cj), Np) = I [e-^-/^ + ai?(6)5(0, b)e 



(32) 



-Et/T' 



.(33) 



The strength of the thermal partons in the medium, char- 
acterized by C'{Np), is not calculable in our approach 
and has always been determined by fitting the \ow-pT 
hadronic data [2l| . Its dependence on Np has been quan- 
tified in [205 based on older data. We now give a better 
determination of it in light of the newer data that include 
also the dependence [l|. 

It has been an assumption in our approach that the 
dependence of pi{pT,<t>j Np) is borne by R{pT,4',Np) 
alone in Eq. and not by B{pT,Np) at all. While 
that assumption has been tested indirectly by V2{pT,b) 
in Fig. 4, we now confront it directly by recognizing that 
if it is valid, then f{pT,4>,Np) in Eq. which we 

can calculate, must correctly describe the experimental 
dependence of pi{pT, 0, Np)^ apart from an overall nor- 
malization. In other words, C{Np) in Eq. ([5^ must be 
independent of 0, a property that could not have been 
checked prior to any knowledge about the dependence 
of /Oi(pti0, Np). With the data on the latter being now 
given in [l|, we have checked that C{Np) is indeed insen- 
sitive to 0. We exhibit that finding in the following more 
revealing way. 
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FIG. 6: The points are empirical C{Np) described in the 
text for pt~1-21 GeV/c (solid circles) and 1.71 GeV/c (open 
squares). The straight line is a power-law fit. 



If we use the experimental values for pi on the LHS of 
Eq. and the calculated values of f{pT, 4>^ ^p) on the 
RHS, we can determine their averages over (p separately; 
we then take the square root of their ratio and regard 
the result as the empirical C{Np). Strictly speaking, it is 
possible for C{Np) thus extracted to depend on whether 
Pt is 1.21 or 1.71 GeV/c. They are, however, shown to be 
independent of pt in Fig. 6 by points depicted in different 
symbols that essentially overlap. We can fit them very 
well by a simple power law 



1.1 iV^-^" GeV" 



(34) 



as shown by the straight line. We can then use this 
formula in Eq. ([32]) and exhibit the (j) dependence of 
Pi{pt, (j^T Np) in t'wo ways. First, we plot the experimen- 
tal data of Pi{pt, 4>: ^p) vs C'^{Np) as points in Fig. 7 at 
the six values of Np treated as parametric variables and 
for (a) 1.0 < PT < 1.5 GeV/c and (b) 1.5 < pr < 2.0 
GeV/c. In each panel there are three sectors of values 
corresponding to = (0, 15°), (30°, 45°) and (75°, 90°) 
analyzed in [i| , although the legend in the figure approx- 
imates them as integral fractions of tt. The curves are 
the results of our evaluation of f{pT,(t),Np) at the cor- 
responding values of the parameters. The agreement be- 
tween theory and experiment is excellent. 

Another way is to exhibit directly i?^^(pTj 0: ^p) 
vs Np in Fig. 8 and RAAipr, (^t ^p) in Fig. 9. 
The calculated results based on Eq. ([32]) are shown 
by solid, dash-dotted and dashed lines for (j) = 
(0,15°), (30°, 45°) and (75°, 90°), respectively They 
agree with the data essentially all within errors. The 
spread for different values of is therefore to be inter- 
preted as a consequence of the physics of ridges. We 
now see that the flatness of the NMF Raa along the 
normal to the reaction plane, pointed out especially in 
[l|, is not caused by any extraordinary, single piece of 
physics. It is due to the cancellation of the linearly ris- 
ing behavior of R'aa shown by the dashed lines in Fig. 
8 and the decreasing behavior of Np/Nc in the rescaling 
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FIG. 7; Pion inclusive distribution pi{pT, 4>, Np) vs C^{Np) 
at six points of A^'part given in [l[| for three ranges of (j) • 
(0,15°), (30°, 45°) and (75°, 90°). (a) 1.0 < pr < 1.5 GeV/c, 
(b) 1.5 < Pt < 2.0 GeV/c. The lines are calculated values 
from Eqs. ^ and (gSJ. 



factor Z{Np)/Z'{Np) in Eqs. (H) and (P. The question 
is why R\j^{<f) = 7r/2, Np) rises with Np. At low px path 
length is not an important issue. The rising behavior 
is the consequence of a combination of factors that in- 
clude the broadening of the tip of the ellipse, when h is 
decreased, and the increasing likelihood for semihard par- 
tons directed at large (pi to produce ridge particles with 
a narrow cone of angular correlation around 02 ~ 7r/2. 
That likelihood is quantified by S{4>, b), as shown by the 
dashed line in Fig. 3. 
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FIG. 8: Alternative R'aa{Pt, 4>y Np) vs Np for (j) and pr values 
as in the figure caption of Fig. 7. Data points are obtained 
from Pi by rescaling Raa{j>t, 0, Np). The lines are calculated 
results. 

In the above we have considered only pion produc- 
tion at low Pt, where TT recombination dominates. For 
Pt > 2 GeV/c, the contributions from TS and SS re- 
combination must be added, as have been calculated 
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FIG. 9: Standard Raa{pt, (j), Np) vs Np. Data points are 
from [l[ as in Fig. 8. The lines are calculated results. 
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FIG. 10: Inclusive distributions of vr" for all centralities rang- 
ing from 0-10% (top) to 80-90% (bottom) in 10% steps, each 
displaced by a factor of 0.2. The data are from [l|]. 



before for all centralities [22]. But now with the new 
parametrization of C{Np) given in Eq. for TT re- 
combination, we can revisit the single-particle distribu- 
tions at higher pr, averaged over all 4>^ even though the 
TT component is not dominant. The calculated result is 
shown in Fig. 10 for nine bins of centrality. The agree- 
ment with data |2l] is evidently very good. How they 
depend on (j) remains to be investigated. 

With this section we have completed the circuit and 
returned to the original question about the characteris- 
tics of the nuclear modification factor raised in Sec. I. 
By reproducing the data on its 4> dependence at low pT, 
we have demonstrated that all dependencies in ridges, 
elliptic flow and NMF are interconnected. 



VIII. CONCLUSION 

We have given a unified description of all the azimuthal 
dependencies of all observables on 7r° production at low 
Pt in heavy-ion collisions. They are: Raa(4>, -^p); V2{Np) 
and ridge yield Y^{(l)s,Np) as functions of Np, the num- 
ber of participants. The main physics input is that the 
semihard scattering near the surface drives the azimuthal 
anisotropy on the one hand and the production of ridge 
particles on the other with or without trigger. The ge- 
ometrical factor, S{(t),b), that makes precise the bridge 
between the two aspects of the problem follows from 
a study of the correlation between the </> directions of 
the trigger and ridge particles. The main understanding 
achieved in this picture is that the single-particle distri- 
bution dN^j^/pTdpTdcj) has two components: one is the 
(/)-independent bulk B{pT,Np), different from the con- 
ventional bulk that is 0-dependent, and the other is the 
ridge component R{pt, 4>j ^p) that carries all the de- 
pendence. It is the latter that can be tested by correla- 
tion experiment using triggers. 

A number of parameters are involved in this analysis, 
but most of them are adopted from previous studies. The 
correlation width a is an important input, and its value 
was determined in Ref. 9] in connection with the study 
of ridgeology. The values of T and T' are from Refs. 
and [2l]| , and that of Sr is an approximation of what 
was found relevant in Ref. The parametrization of 
C{Np) is an improvement of that found in Ref. [l^ and 
is a description of the dependence on centrality, not on (j>. 
The one parameter that we have adjusted in this paper 
to fit the normalization of V2{pTib) is a, which appears 
in Eq. (PT|) and affects the magnitude of the ridge. The 
basic properties of the (j) dependence that we can calcu- 
late, such as S{4>,b) and Y^{(l)s,b) do not depend on a, 
although Raa{ptt4>^^) depends on it indirectly. 

The concern in this paper has been exclusively about 
\ow-pT pions at Pt < 2 GeV/c. Conventionally, such par- 
ticles are regarded as products of soft processes. How- 
ever, our view is that some of them can be due to semi- 
hard processes that can copiously produce semihard par- 
tons whose footprints are the ridge particles at low px- 
Triggers select a subset of events for the study of corre- 
lations, but without triggers the effect of the semihard 
partons must nevertheless be taken into account in the 
calculation of inclusive distributions at low pT- 

The production of proton has not been considered, 
either in this paper or in the experimental study of 
the (j) dependence of Raa [H- Theoretically, large 
baryon/meson ratio has been a signature of the recom- 
bination model [lol |Tl|. Since the same hadronization 
mechanism applies to both inclusive and ridge produc- 
tion, we expect that Raa for proton will have the same 
characteristics as that for pion. There are some differ- 
ences in the pt dependence because of the difference in 
the recombination functions, but the physics connecting 
Raa and is the same. That is, at a fixed centrality 
S{4>, b) exhibits the range of variation when (f) is varied 
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from to 7r/2, which in turn describes the ranges of vari- 
ation of Y^{(j)s) and Raa{4>)j ^ shown in Figs. 5 and 9, 
respectively. The behavior for proton production should 
be similar. This is a good testing ground for any ex- 
otic model that attempts to explain the so-called "baryon 
anomaly" . 

At intermediate pT above 2 GeV/c the dominance of 
TS recombination changes the physics considered here. 
Ridge is a manifestation of the enhanced thermal par- 
tons, and the peak sitting above the ridge is localized in 
A77 and is more closely associated with the hard parton 
through the intermediary shower partons [ll|. Correla- 
tion in a jet peak is different from that in the ridge, so 
the geometrical effect is also different. Although the data 
on Raa{pt, Np) vs A(f> in [l[ do not show significant 
difference between 1.0 < px < 1-5 and 4.0 < pr < 5.0 
GeV/c sectors, their Np dependencies are quite differ- 
ent. For example, Raa{A4' ~ ""/S) is not constant in Np 
for 4 < < 5 GeV/c, in contrast to the behavior at 
1.0 < pt < 1.5 GeV/c. It is clear that the physics at 
intermediate px is very different and requires a fresh ap- 



proach that is not a simple extrapolation from the study 
described here. However, the TT component at higher 
Pt, though smaller, can be reliably calculated only if the 
formalism developed here is used for the description of 
the (j) dependence. Such a study has been carried out 
and a remarkable scaling behavior has been found in the 
dual dependence on and centrality [13], a property that 
is contained in the PHENIX data |l|. 
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